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abstract 


A  frequency  dependent  detection  and  estimation  procedure 
based  on  the  generalized  likelihood  principle  is  applied  to 
several  models  appropriate  for  vertical  and  horizontal  seis¬ 
mic  arrays.  When  data  is  constructed  so  as  to  conform  to 
the  idealized  model  excellent  agreement  with  theoretical  de¬ 
tection  limits  is  shown.  Multiple  signal  models  are  analyzed 
in  the  frequency  domain  with  the  presence  of  various  signal 
components  indicated  in  terms  of  power  spectral  ratios.  Tables 
of  the  central  and  non-central  F  distribution  are  used  to  de¬ 
termine  false  alarm  and  signal  detection  probabilities. 

Two  models  were  assumed  in  analyzing  data  from  vertical 
arrays.  A  full  model,  assuming  each  time  series  to  be  compo¬ 
sed  of  a  ghosted  signal  and  an  infinite  velocity  noise  compon¬ 
ent,  indicated  that  in  general  the  infinite  velocity  component 
is  not  present  in  the  "signal  plus  noise"  or  the  "noise  only" 
cases.  A  reduced  model  containing  only  a  ghosted  signal  plus 
noise  accounted  for  as  much  power  and  produced  as  good  a  sig¬ 
nal  estimate.  The  signal  estimates  in  general  tended  to  be 
point  by  point  reproductions  of  the  surface  trace  indicating 
that  a  significant  signal  to  noise  enhancement  will  be  diffi¬ 
cult  using  a  simple  processor.  Detection  results  indicate 
that  in  "signal  plus  noise"  cases  with  a  high  signal  to  noise 
ratio  the  F  statistic  exceeds  a  threshold  which  would  occur 
by  chance  much  less  than  one  time  in  2000.  However,  "noise 
alone"  analysis  also  gives  a  significant  (although  much  less 
so)  value  over  all  frequencies  indicating  the  presence  of  P 
waves  in  the  noise.  Further  analysis  on  other  arrays  where 
the  fundamental  Rayleigh  mode  is  known  to  be  prominant  and 
where  there  are  lower  signal  to  noise  ratios  should  estat- . 
lish  the  usefulness  of  the  detector.  Analysis  of  the  statis¬ 
tics  of  the  ambient  P-noise  can  also  be  performed. 


For  the  horizontal  array,  results  using  real  data  indi¬ 
cate  that  the  F  statistic,  can  be  a  more  sensitive  detector 
than  the  phased  sum  alone.  Examples  for  the  "signal  plus 
noise"  case  indicate  that  the  signal  is  detected  with  a  high 
probability.  The  "noise  alone"  analysis  yields  a  low  false 
alarm  rate  in  the  particular  cases  analyzed.  Further  research 
can  es  .ablish  the  effectiveness  of  the  statistic  in  determin¬ 
ing  which  frequencies  propagate  with  the  signal  velocity  as 
well  as  the  level  of  detection  possible  with  errors  in  the 
time  delays  and  sensor  calibrations.  Models  which  combine 
vertical  and  horizontal  arrays  are  yet  to  be  investigated. 


INTRODUCTION 


In  the  analysis  of  seismic  array  data  a  natural  exten¬ 
sion  to  the  process  of  estimating  a  signal  by  beamforming  or 
deghosting  is  the  development  of  optimum  measures  for  detect¬ 
ing  the  presence  of  certain  signals  with  a  given  specified 
probability  for  a  false  alarm.  Classical  procedures  involve 
comparing  the  peak  to  peak  amplitude  of  the  processed  trace 
to  the  root  mean  square  of  the  preceding  noise.  Such  proce¬ 
dures  yield  reasonable  results  in  that  for  high  signal  to 
noise  ratios  the  signal  detection  capabilities  are  often 
reasonably  high.  Unfortunately  a  false  alarm  probability  of 
.01  for  a  50-second  sample  would  not  be  unusual  for  such  a 
detector  so  that  the  number  of  false  alarms  in  events  per 
day  becomes  unacceptably  large.  Therefore  the  use  of  ratios 
of  amplitude  to  root  mean  square  noise  levels  must  be  supple¬ 
mented  by  analyst  judgement. 

It  is  interesting  to  note  that  optimum  frequency  de¬ 
pendent  detectors  can  be  developed  in  terms  of  the  general 
likelihood  ratio  criteria.  The  optimum  detectors  involve 
the  comparison  of  various  residual  power  levels  computed 
under  the  assumption  that  various  competing  signal  models 

are  valid.  In  this  report  the  F  detector  is  developed  for 
* 

signal  and  noise  models  which  are  appropriate  for  analyzing 
short  period  horizontal  and  vertical  arrays.  Using  simula¬ 
ted  data  we  establish  that,  for  predictable  false  alarm 
rates ,  propagating  waves  can  be  detected  as  a  function  of 
frequency  with  a  high  level  of  reliability.  When  the  detc- 
tors  are  applied  to  real  data  it  is  noted  that  the  F  statis¬ 
tic  can  isolate  the  frequencies  moving  at  the  signal  velocity 
more  effectively  than  other  measures  such  as  phased  sum  power 
spectra.  The  methods  are  applied  to  two  events  recorded  on 
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a  vertical  array  and  to  two  noise  samples  from  the  same  ver¬ 
tical  array.  A  single  event  recorded  at  the  21  sub-array 
centers  of  LASA  is  analyzed  using  the  theoretical  signal 
model  applicable  to  the  horizontal  array. 


THEORETICAL  DISCUSSION 


Most  signal  models  arising  in  seismo logical  applica¬ 
tions  can  be  expressed  in  the  general  form  (see  Shumway  and 
Dean  [9]) 

Vt}  =  K  I  Xjm(t-u)Bm(u)  +  nj(t) 

J  m=l  u=-°°  J 

j  =  1,  2,...,  N ,  t  =  0,  +_  1,  2,.... 

The  times  series  Y^(t)  is  the  recording  at  jth  sensor  of  a 
vertical  or  horizontal  array.  The  jth  series  Y..(t)  is 
assumed  to  be  a  discrete  time  stationary  process,  and  is 
composed  of  a  mixture  of  signals  S-^Cu),  02 ^ » • • • (u ) .  The 
coefficient  matrix  of  time  invariant  functions  X^Ct)  is 
given  by  specifying  a  particular  model  for  the  signals 
assumed  to  be  present.  The  error  series  n^(t)  are  assumed 
i j  be  uncorrelated  from  sensor  to  sensor.  The  correlation 
function  may  have  a  general  form  although  it  is  assumed  to 
be  identical  from  sensor  to  sensor. 

In  the  next  sections  we  will  present  examples  of  the 
various  special  cases  of  equation  (1)  which  are  relevant 
to  problems  in  seismological  estimation  and  detection.  The 
approach  utilizes  finite  Fourier  transform  methods  to  de¬ 
velop  the  estimates  for  the  signal  and  to  test  the  various 
signal  models  following  the  approach  described  in  refer¬ 
ences  [8],  [9],  and  [10], 

The  Horizontal  Array 

In  particular  we  may  obtain  a  model  for  a  horizontal 
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array  using  equation  (1)  with  p  =  1,  $1(u)  =  s(u)  and  X.^ 

( “t )  ~  6 

t-Tj  where  6^  is  the  Kronecker  delta  function  defined 
to  be  1  when  t  =  0  and  0  when  t  i  0  and  T, ,  T0...,T  are  the 
time  delays  defined  by  the  velocity  and  azimuth  of  the  propa¬ 
gating  signal.  In  this  case  equation  (1)  becomes 


Y.(t)  =  s(t-T.)  +  n.(t)  (2) 


and  we  have  the  signal  model  for  the  horizontal  array. 

Thus  the  problem  reduces  to  the  detection  of  a  determin¬ 
istic  signal  embedded  in  random  noise.  If  we  assume  that  the 
signal  waveform  is  known  and  that  n.(t)  is  a  continuous  para- 
meter  noise  process  possessing  a  Karhunen-Loeve  expansion  in 
terms  of  an  infinite  set  of  uncorrelated  random  variables,  an 
elegant  likelihood  solution  due  originally  to  Grenander  [3] 
is  possible.  However,  the  integral  equations  which  must  be 
solved  in  order  to  generate  the  random  variables  of  interest 
are  difficult  except  for  a  few  simple  cases.  The  assumption 

that  n(t)  is  a  normal  white  noise  process,  for  example,  leads 

r  7 1 

to  the  well  known  matched  filter  (for  one  example  see  Selin  ). 

r  7 1 

Another  possible  approach  can  be  made  via  the  sampling  theorem 
which  also  leads  to  a  matched  filter  detector  for  the  normal 
white  noise  case.  The  difficulties  arising  in  the  solution 
of  the  integral  equations  obtained  via  the  likelihood  approach 
to  the  testing  of  the  hypothesis  s(t)  =  0  in  (2)  against  a 
specified  alternative  have  led  investigators  to  the  use  of 
alternative  methods  for  generating  observable  test  statistics. 

For  example,  using  Fourier  transform  methods  one  can  generate 
random  variables  which  are  uncorrelated  as  the  lengths  of  the 
time  series  tend  to  infinity.  In  addition,  the  sample  vari¬ 
ances  of  such  complex  random  variables  have  expected  values 
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equal  to  the  power  spectrum  if  it  can  be  assumed  that  the 
time  series  is  wide  sense  stationary. 

In  this  section  we  shall  consider  a  general  asymptotic 
solution  to  the  detection  problem  for  a  signal  model  given 
in  the  form  of  equation  (2).  In  particular  we  consider  a 
frequency  dependent  livelihood  ratio  test  of  the  hypothesis 
that  s(t)  =  0  against  the  general  alternative  that  s(t)  t  0. 
The  asymptotic  distribution  of  the  test  criterion  is  given^^ 
in  terms  of  the  classical  F  statistic  enabling  us  to  set  crit¬ 
ical  threshold  values  which  specify  desired  false  alarm  rates. 
Then,  assuming  that  we  have  some  basic  information  on  the 
types  of  alternatives  to  be  tested,  a  procedure  is  presented 
for  computing  the  signal  detection  probability  as  a  function 
of  the  number  of  time  series  and  the  signal  power  in  some 
frequency  band  of  interest.  A  number  of  computational  ex¬ 
amples  verify  that  the  test  statistic  gives  results  which 
are  consistent  with  the  predicted  false  alarm  and  missed 
signal  rates. 

In  particular,  it  has  been  shown^10^  that  the  maximum 
likelihood  estimator  for  the  signal  is  the  phased  sum  of 
the  traces  and  that  the  likelihood  ratio  test  of  the  hypo¬ 
thesis  s(t)  =  0  in  the  frequency  domain  can  be  derived  from 
the  finite  Fourier  transforms  of  the  original  series  Y^Ct) 
and  the  likelihood  estimates  for  the  signal 

1  N 

Y(t)  =  A  l  Y.Ct-T.).  (3) 

j  =  l  J  J 

The  finite  Fourier  transform  is  defined 

T-l 

Y.(wn)  =  i-  T  Y.(t)  u  = 

3  t|0  3  n 


as  usual  by 

— -,  n=  0 ,  1,...,T-1  (4) 
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with  the  transforms  of  Y(t)  and  p^t)  given  respectively  by 
Y(wn)  and  S(to^).  In  [9]  it  is  shown  that  an  analysis  of 
power  given  in  the  form  of  Table  I  below  yields  a  likelihood 
ratio  test  for  the  presence  of  the  signal  with  the  components 
of  power  following  chi-square  distributions. 


TABLE  I 

Analysis  of  Power  at  Frequency 
a)  for  Horizontal  Array 


SOURCE 

POWER 

DEGREES  OF  FREEDOM 

Due  to  Signal 

N 

_  2 
|YOn)| 

2 

N 

2 

Deviation 

I 

|Vj(<on)-Y(<on)| 

L 

2  (N-  1) 

Total 


2N 


W'  refer  to  the  analysis  of  the  variance  involving  the  com¬ 
plex  Hermitian  forms  in  Table  1  as  an  analysis  of  power.  How 
it  can  be  seen  from  Table  I  that  the  power  due  to  the  signal 
is  N  times  the  power  spectrum  of  the  phased  sum  of  the  observed 
time  series  whereas  the  total  power  is  the  sum  of  the  power 
spectra  of  the  original  series.  Since  the  test  statistic  is 
determined  by  a  monotonic  function  of  the  ratio  of  the  power 
due  to  the  signal  to  the  deviation  power,  a  convenient  ratio 
of  two  chi-square  variates  divided  by  their  degress  of  free¬ 
dom  is 


which  gives  a  likelihood  ratio  critical  region.  The  perfor¬ 
mance  of  the  likelihood  detector  above  is  dependent  upon  the 
distributional  properties  of  the  F  statistic  in  equation  (5) 
under  the  null  and  alternative  hypothesis.  If  the  hypothesis 

s(w  )  =  0  is  true  then  (5)  follows  a  central  F  distribution  so 
n 

that  a  specified  false  alarm  rate  (Type  I  error)  a  determines 

an  upper  threshold  value  for  F,  say  F^  ^(2,  2N-2).  For  a 

given  false  alarm  rate  «  the  signal  detection  probability  Pd 

is  determined  by  relative  frequency  with  which  the  observed  F 

will  exceed  F,  (2,  2N-2)  where  the  observed  F  follows  a  non- 

1-a 

central  F  distribution  with  the  noncentrality  parameter  pro¬ 
portional  to 

0  N  I  s  (u)  )  |  2 

62  =  ? -  (6) 

(oMu  )/2 
n 

when  the  signal  is  present.  Mere  o  (wn)  is  the  value  of  the 
theoretical  noise  spectrum  at  frequency  If  the  power  of 

the  signal  in  the  frequency  band  of  interest  is  known  as  well 
as  the  power  of  the  noise,  we  can  compute  a  value  for  the  non¬ 
centrality  parameter  given  in  (6).  Then,  we  determine,  for  a 
given  a,  the  signal  detection  probability  using  the  Fox  charts 
for  the  power  of  the  F  test  (Scheffe  [6],  pp  446-455).  These 
charts  yield  signal  detection  probabilities  for  the  false 
alarm  probabilities  0.01  and  .05  with  various  values  of  the 
noncentrality  parameter.  Since  the  noncentrality  parameter 
is  determined  in  terms  of  the  signal  to  noise  ratio,  it 
follows  that  receiver  operating  characteristic  curves  can  be 
computed  for  the  F  ratio  detector. 

In  theory,  this  procedure  guarantees  a  higher  signal 


detection  probability  than  that  which  would  be  obtained  by 
comparing  the  peak  signal  amplitude  against  a  long  term 
average  noise  level.  In  practice  one  would  prefer  false 
alarm  probabilities  lower  than  .nl,  say  .001,  which  for  a 
60  second  data  sample  would  yie  d  about  1.4  false  alarms 
per  day.  Hence,  for  applications,  the  non  central  F  tables 
would  have  to  be  extended  to  cover  small  false  alarm  prob¬ 
abilities  . 


The  vertical  Array 

For  the  vertical  array  we  assume  a  string  of  vertically 
oriented  sensors  monitoring  a  signal  which  progress  up  the 
array,  is  reflected  off  the  top,  and  then  moves  back  down 
the  array.  In  the  case  of  Rayleigh  and  Love  waves  it  is 
reasonable  to  assume  the  presence  of  a  coherent  noise  trace 
propagating  up  the  array  with  an  infinite  velocity.  If  we 
assume  that  the  relative  amplitude  of  the  noise  varies  in  a 
predictable  fashion  a  general  model  of  the  vertical  array 
becomes 

Y  .  ( t )  =  j  (s ( t-T . )  +  sCt+T.)}  +  a.nCt)  +  n.(t) 

12  j  t-  J  J  J 

j  =  1,  2,  .  .  .N 

where  T, ,  T„,...T  are  the  same  time  delays  characterizing 
the  signal  at  the  various  sensor  locations.  Hence,  each 
time  series  is  composed  of  a  fixed  noise  function  n(t)  and 
a  signal  s(t)  which  appears  as  a  ghost,  first  at  t-Tj ,  and 
then  at  t  +  T..  Equation  (7)  can  be  obtained  as  a  special 
case  of  equation  (1)  by  ^aking  P  =  2, 
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(7) 


81 C  t )  = 


s(t),  $2(t)  =  n(t),  X^(t)  = 


^  ^t-T . + 

3 


't  +  T.  ) 
3 


and  X.«(t)  =  a.  6  where  6,  is  the  Kronecker  delta  as  before. 
3  2  3  t  t 

The  a^  are  constants  which  may  be  thought  of  as  the  amplitude 
of  a  Rayleigh  mode  at  a  particular  frequency. 


In  the  multiple  signal  of  signal-coherent  noise  model 
specified  by  (7)  a  reasonable  sequence  of  tests  of  hypotheses 
might  be  given  by 


H,:  n(t)  =  0 

s(t)  =  0  given  n(t)  i  0 


(8) 


The  above  procedure  tests  first  to  see  if  there  is  an  infi¬ 
nite  velocity  coherent  noise  trace  present.  If  the  hypothesis 
n(t)  =  0  is  rejected  (the  usual  case)  the  processor  tests  the 
hypothesis  that  with  coherent  noise  present  there  is  no  sig¬ 
nal  propagating  up  the  array  with  a  finite  velocity  described 
by  the  time  delays  given  in  (8).  In  [8]  it  is  shown  that  a 
sequence  of  likelihood  ratio  tests  in  the  form  (8)  yields 
the  analysis  of  power  table  given  below  where  the  exact  ex¬ 
pressions  for  the  power  components  are  more  complicated. 
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TABLE  II 


Analysis  of  Power  at  Frequency  u>n  for  a  vertical 

array 


Source 


Power  Degrees  of 

Symbol  Freedom 


Due  to  Coherent  Noise 

Difference  (a)  _ _ 

Due  to  Signal  and  Coherent 
Noise 

Due  to  Incoherent  Noise 


P 

nn 

Diff (a) 


P 

P 


ss  +  nn 
ee 


2 

2 


4 

2 (N- 2) 


Total  Power 


p  2N 

ss+nn+ee 


The  statistic  for  testing  for  coherent  noise  is 

Pnn 

F(  2 »  2H-2)  =  (N-l)  (Dif^aV  ♦  P^T 

which  has  the  F  distribution  with  2  and  2N-2  degrees  of 
freedom.  A  statistic  for  testing  for  the  signal  is  given 

by 


F( 2 ,  2U-4) 


Diff (a) 

TTe 


(N-2 ) 


(10) 


which  has  the  F  distribution  with  ^  and  *.  degree^ 
freedom. 
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For  this  case  the  false  alarm  rate  may  be  specified  as  a 
by  using  F;  (2,2(N-2))  as  a  threshold  value  to  be  exceeded 
only  100a  percent  of  the  time  when  there  is  no  signal  present. 

The  above  detection  procedure  yields  as  residual  results 
asymptotically  best  linear  unbiased  estimates  for  the  signal 
s(t)  and  the  coherent  noise  n(t).  The  estimates  are  also 
asymptotic  maximum  likelihood  estimates.  By  asymptotic  we 
mean  that  the  properties  hold  with  T  large  where  T  is  the 
number  of  data  points. 

The  next  section  presents  several  simulated  examples  de¬ 
signed  to  illustrate  the  estimation  procedure  and  to  verify 
the  performance  of  the  F  statistic  as  an  optimum  detector. 

SIMULATED  EXAMPLES 


In  order  to  determine  whether  or  not  the  finite  versions 
of  the  previously  discussed  asymptotically  optimum  procedures 
give  reasonable  results,  a  number  of  test  cases  were  run 
using  simulated  horizontal  and  vertical  array  data.  The 
results  are  presented  in  the  following  two  sections. 

A  Simulated  Horizontal  Array 

In  order  to  illustrate  the  preceding  sections  it  is 
necessary  to  generate  rets  of  time  series  conforming  to  the 
signal  plus  noise  model  specified  in  equation  (3).  A  128 
point  test  signal  of  the  general  form 

s(t)  =  at  e ”at  sin  8t,  t  =  0,1, ...,127  (11) 

was  generated  where,  as  a  frame  of  reference,  it  is  assumed 
that  the  data  are  obtained  at  a  sampling  rate  of  20  points 


I 


per  second.  Hence  the  frequency  range  covered  runs  from  0 
to  10  cycles  per  second  (cps).  Figure  1  shows  a  128  point 
two  cycle  per  second  signal  generated  according  to  (11). 

This  signal  is  to  be  buried  at  aiffering  magnitudes  in  a 
1024  point  multivariate  stationary  noise  process.  Two  types 
of  noise  processes  were  used  in  the  simulation,  white  Gaussian 
(normal,  independent,  zero  mean)  noise  with  unit  power  density 
at  all  frequencies  and  a  smoothed  white  noise  process  with 
smoothing  coefficients  .10,  .20,  .-10,  .20,  .50  and  -.10.  An 
estimate  for  the  power  spectrum  of  the  smoothed  series  is  shown 
as  the  total  power  component  in  Figure  7.  Figure  2  shows  a 
typical  portion  of  a  1024  point  smoothed  white  noise  series 
containing  a  test  signal  whose  maximum  amplitude  is  0.5.  The 
double  arrow  specifies  the  entry  point  of  the  signal.  For  the 
first  example  ten  series  were  generated  with  each  of  the  two 
types  of  noise  structures  and  with  maximum  signal  amplitudes 
of  0.0,  0.5,  and  1.0  respectively.  Figure  3  shows  the  mear. 
value  of  ten  processes  of  the  same  type  as  that  shown  in  Fig¬ 
ure  2  and  we  see  that  there  is  in  the  averaged  trace  some 
indication  of  a  signal.  A  detector  based  on  signal  ampli¬ 
tude  alone  would  obviously  produce  a  number  of  false  alarms. 
Figures  4,  5,  and  6  summarize  the  analysis  of  power  and  F 
tests  for  the  maximum  signal  amplitude  levels  of  0.0,  0.5, 
and  1.0  respectively  when  the  signals  are  buried  in  white 
Gaussian  noise  with  a  unit  spectrum.  For  practical  applica¬ 
tions  it  is  desirable  to  smooth  the  spectral  estimates  over 
bands  of  approximately  constant  spectral  power.  This  in¬ 
creases  '  he  degrees  of  freedom  for  the  estimates  and  reduces 
the  voluim5  of  data  for  display.  In  effect  this  modifies  Table 

I  in  that  if  the  power  component  for  a  single  n  is  distributed 
2  2 

as  g  (un)  X  9Kr  then  the  average  of  k  identical  independent 

5  1  2  2 

spectral  estimates  would  be  a  (wn)X  2Nk  so  that  entry  should 

2k 
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be  divided  by  N  and  compared  with 


X22Nk 

2Nk 


to  verify  the  theoret- 


2 

ical  distribution.  Therefore  for  a  (u^)  =  1  as  in  Figures  4, 
5,  6a  multiplier  of  2  has  already  been  divided  out  in  the 
smoothing  procedure  so  that  the  observed  chi-square  variables 
in  the  figures  appear  to  have  only  have  as  many  (2,8  and  10) 
degrees  of  freedom  as  would  be  indicated  from  Table  I.  The 
number  of  time  series,  N,  is  ten  for  this  particular  example. 
In  this  set  of  figures  33  estimates  of  the  power  spectrum  are 
made  and  each  estimate  is  an  average  of  17  adjacent  values 
yielding  2  X  17  =  34  degrees  of  freedom  for  the  signal  power 
and  2  x  ^10  -  1)  x  17  =  306  degrees  of  freedom  for  the  error 
power.  If  the  false  alarm  rate  is  set  at  a =  .05,  we  use 
F  q[.  (  30, 300)  =  1.  46  as  the  critical  threshold.  As  the  ampli¬ 


tude  of  the  signal  is  increased  we  may  note  the  improving 
capability  of  the  detector.  When  the  analysis  is  repeated 
using  smoothed  white  noise  (Figures  7,8  and  9)  we  observe 
similar  improvements  in  the  detector  performance.  Table  III 
below  summarizes  some  parameters  of  interest  for  this  first 
example  involving  ten  1024  point  time  series. 


-13- 


TABLE  III 

Parameters  of  Interest  for  10  Time  Series 


Each 

Containing 

a  128  Point 

Signal 

Identification 

Signal 

Maximum . 

Noise 

Non 

m 

Amplitude 

Signal  Power 

Power 

Centrality 

}’:  s’:  i't 

0 

(2.0  cps) 

(2cps ) 

Par ameter 

a  1-8 

Figure  5 

White  Noise 

0.5 

.086 

1.00 

.  9 

.05  .8 

Figure  6 

White  Noise 

1.0 

.  343 

1.00 

1.9 

.01  .99 

Figure  8 

Smoothed 

White  Noise 

0.5 

.086 

.55 

1.3 

.01  .9 

Figure  9 

Smoothed 

White  Noise 

1.0 

.  343 

.  55 

2.5  < 

.01  >.99 

*a  =  Specified  False  alarm  probability 

**1  -  6  =  Pa  =  Theoretical  signal  detection  probability  from 

the  Fox  tables  for  the  power  of  the  F  test. 

The  maximum  signal  power  (in  equation  (6))  is  evaluated  by 
computing  the  average  power  over  the  same  interval  as  the 
spectral  estimates  on  the  basis  of  1024  points.  The  theo¬ 
retical  detection  probability  P^  is  computed  using  the  Fox 
tables  [6  ].  It  is  apparent  from  Figures  4-9  that,  for  a 
variety  of  signals,  a  frequency  by  frequency  detector  based 
on  spectral  power  computations  and  the  F  statistic  gives  re¬ 
sults  consistent  with  theory. 

However,  if  the  procedure  is  to  be  applicable  to  a  broad 
range  of  problems  the  theoretical  false  alarm  and  signal  detec¬ 
tion  probabilities  ought  to  match  the  results  of  simulation 
experiments.  In  order  to  examine  the  performance  measures  over 
a  larger  number  of  trials  and  for  a  limited  number  of  degrees 
of  freedom  a  number  of  simulated  runs  involving  only  three  time 


series  were  made  over  the  parameter  ranges  specified  in  Table 
IV  below.  In  this  case  55  estimates  of  the  power  spectrum 
are  generated  where  each  spectral  estimate  is  an  average  of 
9  adjacent  estimates  yielding  18  degrees  of  freedom  for  the 
signal  power  and  18  x  (3-1)  =  36  degrees  of  freedom  for  the 
error  spectrum.  The  noise  power  at  2  cps  is  .55. 


TABLE  IV 

Theoretical  Parameters  for  Simulation 
(Two  Cycle  per  Second  Signal  Buried  in  Smoothed  White  Noise,  N  =  3) 


Identifi¬ 

cation 

Signal 

Amplitude 

Maximum  Signal 
Power  (2cps) 

Non  Centrality 
Parameter 
(2  cps) 

a 

1-  8 

A 

.145 

.9 

.5 

B 

.327 

1.3 

.7 

C 

582 

1.8 

>  .  9 

A  number  of  repetitions  cf  each  of  experiments,  A,  B,  and  C  in 
Table  IV  above  yielded  the  sample  performance  characteristics 
shown  in  Table  V.  Sample  false  alarm  probabilities  were  es¬ 
timated  by  counting  the  false  alarms  in  the  frequency  bands 
outside  the  2  cps  range. 


TABLE  V 

Sample  Error  Rates  and  Signal  Detection  Probabilities 

For  Simulated  Data 


Identifi¬ 

cation 

Number  of 
Signals 

Number  of 
Signals 
Detected 

Sample  False 
Alarm  Rate 

Sample  Signal 
Detection 
Probability 

A 

34 

15 

.04 

.44 

B 

38 

29 

.01 

• 

cn 

C 

36 

35 

.01 

.97 
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Again  the  sample  results  compare  favorably  with  theoretical 
values  predicted  from  Table  IV. 

A  practical  procedure  then  is  suggested  which  is  similar 
to  the  usual  analysis  of  variance  approach.  First  one  decides 
on  the  frequency  resolution  desired  for  the  analysis  of  power. 
This  determines,  with  a  fixed  data  length,  the  number  of  spec¬ 
tral  estimates  and  the  number  of  degrees  of  freedom  for  each 
spectral  estimate.  As  an  alternative  hypothesis  one  then 
chooses  a  class  of  signals  of  length  L  say  with  a  certain 
waveform  which  are  to  be  detected  with  certain  probability  for 
a  given  false  alarm  rate.  The  noncentrality  parameter  for  the 
F  detector  as  a  function  of  frequency  is  determined  by  appending 
T  -  L  +  1  zeros  to  the  signal  and  then  computing  the  smoothed 
spectrum  in  the  same  fashion  as  for  the  smoothed  estimates  in 
the  analysis  of  power.  For  a  given  false  alarm  rate  the  number 
of  data  series  (sensors)  required  to  detect  the  signal  with  a 
given  probability  can  be  estimated.  If  the  number  of  time  ser¬ 
ies  is  fixed  then  degrees  of  freedom  and  hence  the  signal  de¬ 
tection  level  can  be  adjusted  by  increasing  the  length  of  the 
data  series  or  by  cutting  down  the  number  of  spectral  estimates 
involved  in  each  smoothed  estimaxe. 

A  Simulated  Vertical  Array 

Ten  time  series  containing  1024  points  each  were  genera¬ 
ted  according  to  the  model  specified  by  (7).  In  order  to  pro¬ 
vide  a  check  on  the  probability  distributions  involved  in  the 
analysis  of  power  the  errors  n^(t)  were  taken  as  normal  indepen¬ 
dent  random  variables  with  mean  zero  and  unit  variance.  Under 
this  convention  one  would  expect  the  spectra  of  the  error  traces 
to  be  unity  for  all  frequencies.  A  fixed  coherent  noise  func¬ 
tion  n(t)  to  be  added  to  each  trace  was  generated  by  taking  three 
point  running  averages  of  zero  mean  independent  normal  random 
variables.  This  generates  a  series  with  a  predictable  power 
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spectrum.  A  portion  of  the  1024  point  fixed  function  n(t) 
is  shown  in  Figure  14  and  its  theoretical  spectrum  is  given 
in  Figure  16.  An  arbitrary  sampling  rate  of  20  points  per 
second  is  assumed  which  yields  a  ten  cycle  per  second  folding 
frequency  corresponding  to  the  point  tt  .  The  signal  s(t)  is 
generated  as  a  relatively  short  (128  point)  exponentially  de¬ 
caying  sine  wave  oscillating  at  three  cycles  per  second,  and 
is  shown  in  Figure  12.  Figure  10  shows  a  240  point  portion 
of  one  of  the  ten  time  series  simulated  according  to  the  model 
implied  by  (7)  with  a^=a2= . . . =a1Q=l.  The  maximum  amplitude 
of  the i signal  function  is  4.0  while  the  coherent  noise  and 
random  error  have  unit  variance  each.  The  arrows  signal  the 
two  entry  points  for  s(t)  and  are  respectively  23  points  to  the 
left  and  right  of  the  center  point. 

Now  for  the  vertical  array  model  (7)  with  a^=l  we  may 

write  X.  (t)  =  ^[6  „  +  6  „  ]  where  6.  is  1  when  tt  °  and  0 
1 1  t- 1 .  t+  r .  t 

]  D 

otherwise.  Then  we  may  write  the  Fourier  transform  of  X^(t)  as 


:-jl(a>)  =  |  ^t-T.  +  6t+T.J 


-ioit  m 

e  =  cos  coT. 

: 


(12) 


with  X.,;(io)  =  1  so  that  the  solutions  for  the  best  linear  un- 
biased  filters  in  frequency  may  be  written  using  the  usual 
matrix  equation  [9]. 


H(co)  =  (X*X)' 


X*(to) 


(13) 


for  the  optimum  filters  at  each  frequency.  This  determines 


H,  ,(<*))  =  x  [N  cosojT,  -  b(w)] 
1K  I(Jj)  K 


(14) 


H2k(0 


1  [A(o>)  -  B(oj)  cos  0)T.  ] 

Alu)  k 


(15) 
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with 


N 


N 


A(w)  =  l  cos  wT^,  B(co)  =  l  cos  u)Tv ,  A(co)  =  NA(w)  - 


k=  1 


k=  1 


k» 


B‘  (w) 


(16) 


•/here  k=l,...,N.  In  theory  the  frequency  functions  (14)  - 
(16)  could  be  multiplied  and  summed  frequency  by  frequency  with 
the  transforms  of  the  original  time  series  in  the  frequency  do¬ 
main  version  of  a  matrix  convolution.  In  practice,  we  avoid 
end  effects  and  simulate  the  infinite  two  sided  nature  of  the 
convolution  by  constructing  the  time  functions  hj^(t)  over  a 
shorter  time  interval  than  the  interval  over*  which  the  original 
time  series  y^(t)  are  known.  The  interval  of  length  L  say 
must  be  long  enough  to  cover  the  largest  time  delays  specified 
by  the  linear  model  and  short  enough  to  avoid  the  end  effects 
inherent  in  finite  convolutions.  Then  we  may  compute  the  fre¬ 
quency  functions  at  each  of  the  frequencies  ok  =  27T./L,  j  =  0, 

1 , .  .  .  , L-l .  The  complication  at£(Q)  =  0  is  eliminated  by 

assigning  H.,  (0)  =  0  so  that  the  filters  do  not  pass  the  zero 
3  K 

frequency.  The  finite  Fourier  transform  applied  to  the  fre¬ 
quency  functions  (14)  -  (16)  yielded  the  time  functions  shown 
in  Figure  11  for  N  -  10  and  L  *  128  points.  The  filter  for 
estimating  the  mean  effect  contains  the  expected  peak  at  time 
zero  for  all  h^(t).  The  signal  estimation,  or  s(t)  filters, 
contain  the  peaks  at  the  time  delays  at  v/hich  the  signal  appears. 
For  the  tenth  series  the  time  delay  was  23  points.  The  negative 
small  peaks  match  the  signal  time  delays  at  the  other  levels. 
Hence  we  have  a  matrix  of  filters  reminiscent  of  the  adjusted 
row  and  column  totals  which  estimate  effects  in  the  classical 
analysis  of  variance.  Now  the  128  point  filters  may  be  extended 
by  adding  zeros  and  convolved  with  the  multivariate  time  series 
(Yj(t),  k  =  1,...,10,  t=l, ..., 1024) using  the  usual  finite  Fourier 
■transform  approach. 
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The  true  and  estimated  signals  are  shown  in  Figure  12 
and  Figure  13.  Portions  of  the  true  and  estimated  noise 
functions  are  shown  in  Figure  14  and  Figure  15  respectively. 

The  portion  of  the  estimated  signal  shown  is  a  segment  which 
was  originally  confounded  with  the  noise  function  on  the  original 
traces  and  it  can  be  seen  that  both  estimates  are  remarkable 
repr  ductions  of  their  known  signatures.  It  is  of  interest  to 
point  out  that  the  method  of  estimation  is  independent  of  any 
postulated  waveform  for  either  the  signal  or  the  noise. 

For  the  linear  model  presented  here  we  are  interested  in 
the  following  sequence  of  tests, 

H^:  n(t)  =  0 

H^:  s(t)  =  0  given  n(t)  i  0 

This  is  seen  to  be  an  example  of  a  sequence  of  tests  of  the 
kind  given  in  Table  II.  As  a  practical  matter  the  power  due 
to  n(t)  is  seen  to  be  the  power  spectrum  of  the  predicted 
trace  Y(t)  under  the  assumption  that  s(t)  =  0. 

Thus  the  power  due  to  the  mean  is  N|Y(w)|  which  is  com¬ 
puted  from  the  mean  value  of  the  time  functions  by  the  usual 
finite  Fourier  transform  methods.  Similarly  it  can  be  seen 
that  the  power  due  to  fitting  both  nCt)  and  s(t)  quoted  as 
P  .  in  Table  III  is  just  the  sum  of  the  power  spectra  of 
the  predicted  time  series  when  neither  n(t)  or  s(t)  are  zero. 
This  term  could  be  computed  directly  in  the  frequency  domain 
but  is  more  easily  computed  by  constructing  the  predicted 
series  at  each  level  under  the  assumption  that  neither  s(t) 
nor  n(t)  are  zero  and  then  summing  the  power  spectra  of  these 
predicted  series  computed  by  the  usual  finite  Fourier  transform 
methods.  The  total  power  is  obtained  by  summation  of  the  power 
spectra  of  each  of  the  original  time  series.  Other  terms  are 
obtained  by  differencing. 
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Figures  17,  18  and  19  summarize  the  analysis  of  power 
results.  In  Figure  1 7  we  observe  the  three  successive  power 
components  as  a  function  of  frequency.  In  general  the  com¬ 
putational  procedure  for  a  power  component  here  was  to  compute 
31  point  running  averages  of  the  squared  finite  Fourier  trans¬ 
forms  of  the  original  1024  point  time  series.  This  yielded 
33  power  spectral  estimates  with  approximately  60  degrees  of 
freedom  each.  Adjacent  estimates  are  slightly  correlated.  In 
Figure  17  the  power  due  to  the  noise  plus  signal  function 
shown  as  the  solid  line  contains  the  expected  peak  at  3  cycles 
per  second,  the  frequency  of  the  signal  function  s(t).  The 
difference  between  the  total  power  and  the  power  due  to  s(t) 
is  uniform  at  all  frequencies  as  is  expected  since  the  input 
error  trace  has  a  constant  power  spectrum.  Figure  18  shows 
the  F  statistic  for  testing  the  noise  function  computed  from 
Equation  (10).  The  degrees  of  freedom  are  then  because  of  the 
smoothing  approximately  31x2xl«60  and  31x2  (10-2  )*s500  .  Since 
the  error  spectrum  is  again  uniform  the  F  statistic  as  a  func¬ 
tion  of  frequency  approximates  the  spectrum  given  in  Figure  16. 
Obviously  coherent  noise  is  present.  Figure  19  shows  the  F 
statistic  computed  using  Equation  (10)  for  testing  the  signal 
function  and  we  see  the  expected  highly  significant  peak  at 
the  signal  frequency  of  3  cycles  per  second.  At  other  frequen¬ 
cies  the  statistic  stays  below  its  critical  value.  Thus  when 
the  signal  is  present  even  as  a  128  point  function  buried  in 
a  1024  point  series  the  F  ratio  statistic  is  a  highly  sensitive 
detector . 

In  order  to  examine  che  analysis  of  power  when  the  signal 
was  not  present  the  entire  procedure  was  repeated  with  data 
simulated  according  to  the  signal  free  model.  Again  the  error 
power  has  a  uniform  spectrum  with  the  appropriate  degrees  of 
freedom,  and  the  F  statistic  for  testing  the  significance  of 
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the  signal  function  s(t)  conforms  well  to  its  null  F  distri¬ 
bution  with  60  to  500  degrees  of  freedom.  Of  course,  the  F 
statistic  for  testing  the  noise  'function  n(t)  still  remains 
significant.  Finally  the  analysis  of  power  was  repeated  on 
simulated  data  containing  only  the  error  process  e.(t).  This 
time  neither  of  the  F  statistics  are  more  than  occasionally 
significant  and  they  both  conformed  approximately  to  the  null 
distributions. 

SEISMIC  ARRAY  DATA 


In  the  previous  sections  we  have  established  that  for 
data  conforming  to  one  of  the  two  basic  models,  the  detection 
and  estimation  procedure  presented  in  Section  2  is  effective. 
The  excellent  agreement  between  the  theory  and  the  results  of 
the  simulation  confirms  the  correctness  of  the  computer  pro¬ 
grams  and  the  validity  of  the  asymptotic  approximations  used. 
The  next  logical  step  is  to  examine  the  performance  of  the  op¬ 
timum  detector  on  a  reasonable  test  set  of  seismic  data.  For 
this  purpose  a  horizontal  array  composed  of  21  subarray  centers 
at  LASA  and  a  representative  short  period  vertical  array  of  re¬ 
cording  instruments  including  the  surface  was  chosen. 

Horizontal  Array 

A  representative  event  at  the  21  subarray  centers  for 
LASA  is  shown  for  several  sensors  in  Figure  20.  An  analysis 
of  power  (see  Table  I)  for  this  event  yielded  the  power  com¬ 
ponents  shown  in  Figures  21  and  22.  The  power  spectrum  of  • 
the  beam  formed  sum  shows  significant  power  components  at  .5 
and  1.0  cps.  However  the  F  statistic  indicates  that  the 
signal  model  implied  by  the  beam  is  only  significant  at  . 5 
cps  as  can  be  deduced  from  Figure  23,  Hence  the  power  at  1 
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cps  is  not  explained  by  the  signal  model  and  may  be  due  either 
to  site  reverberations  or  signal  generated  noise.  The  power 
due  to  the  unexplained  noise  (Figure  22)  also  shows  a  strong 
contribution  at  1  cps.  In  order  to  check  the  performance  of 
the  F  statistic  on  noise  the  analysis  was  repeated  when  the 
signal  was  absent.  The  result  shown  in  Figures  24  -  26  in¬ 
dicate  that  at  a  reasonable  threshold  no  false  alarms  were 
observed  for  this  particular  case. 

The  processor  is  useful  not  only  as  an  improved  automatic 
'  detector  but  can  isolate  frequencies  which  are  not  propagating 
with  the  beam  characteristics.  It  is  foirly  simple  to  extend 
the  results  to  the  case  where  more  than  the  signal  is  present 
as,  for  example,  to  the  model 

Y . (t )  =  s. ( t-T . )  +  s0  (t-D.)  +  e. (t) .  (17) 

The  estimation  of  both  signals  has  been  considered  (see  [9]) 
and  the  detection  theory  can  be  developed  for  testing  whether 
one  or  both  of  the  signals  is  present  at  a  given  frequency. 
Hence  the  procedure  will  discriminate  on  a  frequency  by  fre¬ 
quency  basis  between  signals  from  different  sources  or  detect 
the  presence  of  directional  noise.  The  sensitivity  of  the  F 
statistic  as  a  means  for  resolving  signals  from  different 
directions  is  yet  to  be  investigated. 

Vertical  Array 

Previous  investigations  of  the  vertical  array  using  the 
likelihood  filters  and  the  plane  parallel  layering  Raylei.gh 
noise  model  as  described  in  Dean  [2]  have  not  been  effective. 
It  was  established  that  for  simulated  data  Rayleigh  Wave  re¬ 
jection  by  optimum  filtering  would  be  possible.  This  section 
describes  some  further  experimental  results  using  two  simpli¬ 
fied  models  along  with  the  F  statistic  as  a  diagnostic  model 
indicator. 
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For  the  data  analysis,  a  short  period  vertical  array 
(site  207)  recorded  an  event  from  the  Fiji  Islands  and  an 
event  from  the  east  coast  of  Colombia*  Time  de^av  , 
from  the  records  (see  Equation  (7))  were  .00,  «21,  .-8, 

.36,  .M3,  and  .82  seconds  respectively.  These  are  one  way 
travel  tim  s  read  starting  at  the  surface. 

All  data  are  prepared  by  a  bandpass  filter  operating 
over  the  band  (.6-2.0)  cycles  per  second.  Figure  28  shows 
the  basic  data  used  as  input  to  the  optimum  processor  pro¬ 
gram.  A  visual  examination  of  the  various  recorded  outputs 
indicates  that  the  noise  preceding  the  signal  might  conform 
reasonably  well  to  the  noise  model  given  in  (7)  and  for 
reference  we  reproduce  the  basic  vertical  array  model  to  be 

considered  as 

Y.(t)  =  %{  s ( t-TA )  +  s(t  +  T.)}  +  a^n(t)  +^(t),  (18) 

j  =  1.....7 

The  noise  preceding  the  signal  appears  to  have  the  same  wave¬ 
form  from  level  to  level  except  for  the  amplification  factor 
a. ,  but  it  is  difficult  to  determine  the  velocity  of  propaga¬ 
tion  which  may  be  close  to  that  of  the  signal.  The  relative 
amplification  factors  were  estimated  by  examining  a  noise 
record  preceding  the  signal.  The  relative  amplitudes  of  the 
apparent  infinite  velocity  noise  were  read  for  each  event  and 
were  not  assumed  to  vary  over  the  frequency  range  of  the  data. 
The  a . 1 s  did  not  in  general  attenuate  with  depth  as  would  be 
expected  if  n(t)  were  fundamental  mode  Rayleigh  noise.  A 
typical  section  of  noise  record  is  shown  in  Figure  30  for 
14007  and  again  we  notice  the  high  coherence  between  the 
sensor  outputs.  If  a  separation  between  the  signal  and  noise 
propagation  modes  is  possible  then  the  analysis  of  power  given 
in  Table  II  should  show  separate  components  due  to  the  signal 
and  coherent  noise. 
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In  order  to  apply  the  technique  described  in  the  previous 
sections  .  to  the  preceding  vertical  array  data  the  filters 
given  in  Equations  (14),  (15),  and  (16)  were  generated  in  fre¬ 
quency  and  the  results  were  transformed  to  the  time  domain. 

One  of  the  resulting  filters  in  frequency  and  time  is  shown  m 
Figure  27  for  both  the  signal  estimator  and  the  noise  estimator. 
The  signal  estimation  filters  in  time  contain  the  expected  peaks 
at  the  positive  and  negative  time  delays  with  a  negative  noise 
suppression  peak  in  the  center. 

When  the  optimum  filters  were  applied  to  ti  <■  ver:ical 
arrays  with  the  signal  present,  the  filtered  outputs  (Figures 
29  and  33)  give  an  almost  identical  reproduction  of  the  signal 
as  it  appears  at  the  surface.  Since  the  noise  level  compared 
to  the  surface)  is  not  reduced  appreciably  on  the  filtered  out¬ 
put,  there  appears  to  be  significant  noise  propagating  with  the 
signal  velocity.  The  estimated  noise  traces  (Figures  29  and 
33)  indicate  signal  generated  noise  in  the  case  of  Figure  29 
but  not  as  much  signal  generated  noise  in  the  case  of  Figure 
33.  If  the  filters  are  convolved  with  a  set  of  outputs  con¬ 
taining  noise  alone  (Figure  31),  both  filtered  outputs  look 
similar  to  the  apparent  common  noise  trace  as  can  be  observed 
from  Figure  30.  This  indicates  that  significant  energy  in  the 
noise  may  be  propagating  at  both  P  wave  and  Rayleigh  velocities. 

In  order  to  examine  the  apparent  difficulty  of  the  filters 
in  separating  the  signal  velocity  and  noise  velocity  components, 
it  is  interesting  to  look  at  the  detector  results  as  read  from 
the  analysis  of  power  giv>  .  by  using  Table  II.  For  the  analysis 
of  power  computations  the  following  parameters  and  conventions 
are  used.  The  sampling  rate  is  approximately  40  points  per 
second  yielding  a  folding  frequency  of  20  cycles.  Estimates 
are  generated  from  2048  data  points  with  the  spectral  estimates 
composed  by  taking  15  point  running  averages  which  yielded 


15x2  or  30  degrees  of  freedom  per  spectral  estimate.  With  7 
data  series  the  degrees  of  freedom  for  the  error  spectrum  are 
(7-2 )xl5x2  =  150  so  that  the  appropriate  F  statistic  has  30 
and  150  degrees  of  freedom. 


Table  VII  shows  the  analysis  of  power  results  on  the  Fiji 
Island  event  when  the  full  model  given  by  Equation  (18)  is 
applied.  In  this  first  case,  the  signal  is  present  (Figure  28) 
and  the  F  statistic  for  detecting  the  signal  is  appropriately 
significant  over  the  whole  range  at  the  1%  level.  The  maximum 
value  for  the  total  power  component,  computed  by  summing  the 
individual  original  power  spectra,  occurs  at  .8  cps  whereas 
the  maximum  of  the  F  statistic  for  detecting  the  signal  is 
Closer  to  the  usual  signal  frequency  (1.4-1. 6  cps).  We  notice 
here  that  the  F  statistic  for  detecting  coherent  infinite 
velocity  noise  exceeds  the  critical  level  only  at  a  very  low 
frequency  (.4  cps)  implying  that  we  reject  the  hypothesis  that 
such  a  coherent  noise  signature  is  present  except  at  low  fre¬ 
quencies.  It  will  be  noticed  that  in  the  lower  frequencies  the 
power  components  are  not  always  additive  so  that  F  values 
become  large  or  are  negative.  These  F  values  are  included  in 
the  Tables  for  completeness  but  should  not  be  relied  upon  as 
indicators  for  the  presence  of  signal  or  noise  at  low  frequen¬ 
cies.  The  reason  for  the  non  additivity  lies  in  the  restriction 
that  the  filters  not  pass  Love  frequency.  Hence,  for  128  point 
filters,  the  frequency  response  is  incorrect  from  0  to  .16  cps. 
Now  the  spectral  averages  are  computed  by  15  point  running 
averages  which  lead  to  a  bandwidth  of  approximately  13  cps. 

Hence  an  estimate  centered  at  .6  cps  is  about  the  lower  limit 
of  credibility  for  the  power  components  and  F  ratios.  Lower 
frequencies  can  be  examined  by  decimating  the  data  and  taking 
a  longer  time  series.  To  verify  that  the  infinite  velocity 
noise  form  does  not  describe  the  data  adequately,  a  section  of 
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data  immediately  preceding  the  Fiji  Island  signal  was  analyzed 
in  a  "noise-alone  case".  The  analysis  of  power  (Table VII/)  shows 
again  that  the  "signal-like  noise"  dominates  the  noise  which 
might  be  propagating  with  infinite  velocity.  The  F  statistic 
indicates  at  a  high  significance  level  that  a  signal  is  present 
or  that  a  time  trace  with  a  signal  velocity  can  be  isolated. 

And  also,  that  infinite  velocity  noise  is  present.  In  Table  VIII 
a  Colombian  event  is  analyzed  and  the  results  are  similar  to 
the  Fiji  Island  case.  The  F  statistic  for  detecting  the  signal 
is  highly  significant  at  1.6  cps  with  the  noise  F  ratio  only 
occasionally  significant. 

Since  the  noise  propagating  at  infinite  velocity  was 
generally  not  significant,  we  may  consider  a  reduced  version 
of  (18)  given  by 

Y.(t)  =  J5  (s(t-T.)  +  s(t+T.)}  +  nj(t)  (19) 

which  does  not  include  a  coherent  noise  component.  The  results, 
shown  in  Figures  34  -  36  tend  to  confirm  that  the  optimum  fil¬ 
ters  for  the  reduced  model  produce  an  adequate  version  for  the 
estimated  signal.  The  analysis  of  power  using  the  reduced 
model  (19)  is  shown  in  Tab.lv  IX  and  we  see  that  the  assignment 
of  the  coherent  noise  power  lerm  to  the  error  spectrum  has  not 
adversely  affected  the  power  due  to  the  signal.  By  comparing, 
for  example,  the  frequency  band  centered  at  1.4  cps  for  Table 
IX  with  the  corresponding  entries  in  Table  VI  we  see  that  the 
error  power  (39.52  -  30.19  =  9.33)  increases  to  12.59  reflect¬ 
ing  some  of  the  coherent  noise  moving  into  the  incoherent  error 
spectrum. 

The  F  statistic  for  detecting  the  signal  continues  to  be 
highly  significant  in  the  signal  band  (1.4  -  1.6  cps).  In  the 
signal  free  case  (Table  X  )  the  noise  propagating  with  a  signal¬ 
like  velocity  again  gives  a  significant  number  of  false  alarms 
using  the  F  statistic.  However,  the  magnitudes  of  the  observed 
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TABLE  VI 

Analysis  of  Power  for  Fiji  Island  Event  (Full  Model,  Signal  Only) 


t 


TOO 


TABLE  VII 

Analysis  of  Power  for  Fiji  Island  Event  (Full  Model,  Noise  Only) 


F  ratios  are  considerably  lower  than  in  the  case  where  the 
signal  is  present  and  offer  a  reasonable  basis  for  deciding 
whether  a  signal  is  present.  The  fact  that  there  is  very 
little  power  in  the  signal  band  for  Table  X  suggests  that  we 
are  detecting  very  weak  P  arrivals  which  are  apparently  always 
present.  The  analysis  of  power  for  the  Colombian  event  (Table 
XI)  shows  similar  results  and  can  be  compared  with  the  analysis 
given  in  Table  VIII. 

The  above  experimantal  results  tend  to  indicate  that  the 
reason  for  the  unimpressive  filtering  results  is  the  presence 
of  significant  noise  which  propagates  within  the  signal  band 
at  the  same  velocity  as  the  signal.  This  phenomenon  can  be 
examined  by  a  frequency  wave  number  analysis  of  the  noise  alone 
and  the  signal  plus  noise. cases.  Figures  37  and  38  show  the 
frequency  wave  number  spectra  for  data  on  the  Fiji  Island  event 
limited  to  the  frequency  band  1-2  cps  and  we  may  note  again  the 
almost  identical  velocity  characteristics  of  the  noise  and  the 
signal. 

DISCUSSION 


In  this  study  of  optimum  estimation  and  detection  pro¬ 
cedures  for  seismic  arrays  it  was  established  that  when  theoret¬ 
ical  models  postulated  for  vertical  or  horizontal  arrays  hold 
(as  in  the  simulated  examples  in  Section  3)  the  methods  will 
perform  as  advertised.  The  construction  of  intuitively  pleasing 
power  ratios  as  F  detectors  and  frequency  dependent  maximum  like¬ 
lihood  estimation  procedures  present  few  computational  problems. 

An  example  showed  that  the  horizontal  array  detector  can 
aid  in  the  model  specification  problem  by  isolating  signals 
from  identical  or  different  sources.  We  also  were  aDle  to  give 
an  optimum  statistic  for  detecting  the  signal,  expressed  in 
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terms  of  the  individual  spectra  and  the  phased  sum  spectrum. 
Using  this  detector  one  can  set  an  arbitrary  false  alarm  prob 
ability  which  determines  a  signal  detection  proba'  ility  in  terms 
of  the  signal  power  divided  by  the  noise  power.  An  event  re¬ 
corded  at  21  LASA  subarray  centers  gave  reasonable  results. 

A  basic  difficulty  in  the  application  of  the  likelihood 
estimation  procedure  to  the  vertical  array  was  that  F  wave 
noise  was  present  which  interfered  with  the  estimated  signal 
waveform.  This  interference  combined  with  the  degradation  of 
the  signal  waveform  down  the  array  kept  the  signal  to  noise 
ratio  of  the  processed  trace  about  the  same  as  the  surface 
trace.  However,  the  vertical  array  shows  promise  as  a  detector 
even  in  the  presence  of  interfering  P  noise.  The  F  statistic 
shows  some  stability  for  the  "noise  alone"  case  for  the  two 
events  analysed  and  exceeds  a  very  high  threshold  when  the 
signal  is  present.  The  failure  of  F  to  follow  its  exact 
theoretical  signal  free  distribution  when  noise  only  is  pies- 
ent  need  not  indicate  that  the  detection  procedure  will  fail. 
One  may  measure  the  F  associated  with  ambient  P  noise  to  estab¬ 
lish  a  null  non  central  F  distribution.  If  the  F  statistic 
computed  for  the  signal  present  case  exceeds  its  threshold  con¬ 
sistently  evtn  for  low  signal  amplitudes  the  detector  ^an  still 
perform  consistently.  This  is  presently  under  investigation 
as  well  as  are  efforts  to  realistically  model  the  ambient 

noise . 

In  summary,  we  have  shown  that  the  generalised  likelihoc  d 
ratio  approach  can  solve  simultaneously  the  detection  and  esti¬ 
mation  problems  or,  a  frequency  by  frequency  basi  .  A 
general  model  can  be  used  which  includes  many  multiple  signal , 
multiple  noise  propagation  modes  as  special  eases,  ih^  pro- 
*  cedure  is  general  enough  to  be  applicable  to  horizontal  array  , 

vertical  arrays,  or  combinations  of  horizontal  and  vertical 
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arrays.  A  major  application  of  the  generalized  theory  pre¬ 
sented  is  the  testing  of  one  postulated  theoretical  model 
against  another.  In  this  way  the  t na  *ysis  of  power  aids  in 
the  important  task  of  model  specification  using  reasonable 
statistical  testing  procedures. 
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Figure  6.  Analysis  of  power,  Sig.  Amp.  =  1.0,  Noise  Power 
=  1.0,  Number  of  Sensors  =  10. 
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Analysis  of  power,  Sig.  Amp.  =  0.0,  Noise  Power 
=  0.55,  Number  of  Sensors  =  10. 
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Figure  17.  Analysis  of  power. 
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Figure  21.  Estimated  power  due  to  signal  (phased  sum) 
DF  =  6,  signal  is  present. 
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Figure  25.  Estimated  power  due  to  noise  DF  =  120  signal  is 
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Figure  26.  Signal  to  noise  ratio  F  statistic  with  6  and 
120df  signal  is  present. 
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Input  seismic  traces  containing  signal  from 
vertical  array  recording  event  from  Fiji 
Islands  (seismogram  14007). 


Figure  30.  Input  seismic  traces  containing  noise  from 
vertical  array  recording  event  from  Fiji 
Islands  (seismogram  14007). 
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Figure  32.  Input  seismic  traces  containing  signal  from 
vertical  array  recording  an  event  from 
Colombia  (seismogram  14012). 
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Figure  38.  Frequency  wave  number  analysis  with  noise 
present  for  Fiji  Island  event. 
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13  ABSTRACT 

A  frequency  de 
based  on  the  genera 
several  models  appr 
mic  arrays.  When  d 
idealized  model  exc 
limits  is  shown.  M 
frequency  domain  wi 
indicated  in  terms 
central  and  non-cen 
false  alarm  and  sig 
examples  us i ng  real 


pendent  detection  and  estimation  procedure 
lized  likelihood  principle  is  applied  to 
opriate  for  vertical  and  horizontal  seis- 
ata  is  constructed  so  as  to  conform  to  the 
el  lent  agreement  with  theoretical  detection 
ultiple  signal  models  are  analyzed  in  the 
th  the  presence  of  various  signal  components 
of  power  spectral  ratios.  Table  of  the 
itral  F  distribution  are  used  to  determine 
inal  detection  probabilities.  A  number  of 
seismic  data  are  presented. 
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